Low-energy theory and RKKY interaction for interacting quantum wires with Rashba 

spin-orbit coupling 
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We present the effective low-energy theory for interacting ID quantum wires subject to Rashba 
spin-orbit coupling. Under a one-loop renormalization group scheme including all allowed interaction 
processes for not too weak Rashba coupling, we show that electron-electron backscattering is an 
irrelevant perturbation. Therefore no gap arises and electronic transport is described by a modified 
Luttinger liquid theory. As an application of the theory, we discuss the RKKY interaction between 
two magnetic impurities. Interactions are shown to induce a slower power-law decay of the RKKY 
range function than the usual ID noninteracting cos(2fc_Fx)/|a;| law. Moreover, in the noninteracting 
Rashba wire, the spin-orbit coupling causes a twisted (anisotropic) range function with several 
different spatial oscillation periods. In the interacting case, we show that one special oscillation 
period leads to the slowest decay, and therefore dominates the RKKY interaction for large separation. 

PACS numbers: 71.10.Pm, 85.75.-d, 73.63.-b 



I. INTRODUCTION 



Spin transport in one-dimensional (ID) quantum wires 
continues to be a topic of much interest in solid-state and 
nanoscale physics, offering interesting fundamental ques- 
tions as well as technological applications^* Of particular 
interest to this field is the spintronic field effect tran- 
sistor (spin-FET) proposal by Datta and Das^ where 
a gate-tunable Rashba spin-orbit interaction (SOI) of 
strength a allows for a purely electrical manipulation of 
the spin-dependent current. While the Rashba SOI arises 
from a structural inversion asymmetry^*-^ of the two- di- 
mensional electron gas (2DEG) in semiconductor devices 
hosting the quantum wire, additional sources for SOI can 
be present. In particular, for bulk inversion asymmetric 
materials, the Dresselhaus SOI (of strength (3) should 
also be taken into account. By tuning the Rashba SOI 
(via gate voltages) to the special point a = (3, the spin- 
FET was predicted to show a remarkable insensitivity to 
disorder^ see also Ref. 0- On top of these two, additional 
(though generally weaker) contributions may arise from 
the electric confinement fields forming the quantum wire. 
In this paper, we focus on the case of Rashba SOI and 
disregard all other SOI terms. This limit can be realized 
experimentally by applying sufficiently strong backgate 
voltageS ) 8 i 9 i 10 ' 1:L which create a large interfacial electric 
field and hence a significant and tunable Rashba SOI cou- 
pling a. The model studied below may also be relevant to 
ID electron surface states of self-assembled gold chains J* 2 - 

The noninteracting theory of such a "Rashba quantum 
wire" has been discussed in the literature , 13 i 14 ' 15 i 16 ' 17 i 18 
and is summarized in Sec. [IT] below. We here discuss 
electron-electron (e-e) interaction effects in the ID limit, 
where only the lowest (spinful) band is occupied. The 
bandstructure at low energy scales is then characterized 



by two velocities^ 

va.b = i>f(1 ± 8), 8(a) cx a 4 . 



(1) 



These reduce to a single Fermi velocity Vf hi the absence 
of Rashba SOI (8 = for a = 0), but they will be differ- 
ent for a^O, reflecting the broken spin SU (2) invariancc 
in a spin-orbit coupled system. The small-a dependence 
8 cx a 4 follows for the model below and has also been 
reported in Ref. H^- Therefore, the velocity splitting (JTJ 
is typically weak. While a similar velocity splitting also 
happens in a magnetic Zeeman field (without SOI),— the 
underlying physics is different since time-reversal symme- 
try is not broken by SOI. 

The bandstructure of a single-channel quantum wire 
with Rashba SOI should be obtained by taking into ac- 
count at least the lowest two (spinful) subbands, since a 
restriction to the lowest subband alone would eliminate 
spin relaxation. 15 i 22 ' 23 The problem in this truncated 
Hilbcrt space can be readily diagonalizcd, and yields two 
pairs of energy bands. When describing a single-channel 
quantum wire, one then keeps only the lower pair of these 
energy bands. We mention in passing that bandstructure 
effects in the presence of both Rashba SOI and magnetic 
fields have also been studied^^i^i 2 -^ In addition, the 
possibility of a spatial modulation of the Rashba coupling 
was discussed^ but such phenomena will not be further 
considered here. Finally, disorder effects were addressed 
in Refs. H3M 

For ID quantum wires, it is well known that the in- 
clusion of e-c interactions leads to a breakdown of Fermi 
liquid theory, and often implies Luttinger liquid (LL) be- 
havior. This non- Fermi liquid state of matter has a num- 
ber of interesting features, including the phenomenon of 
spin-charge separation.— Motivated mainly by the ques- 
tion of how the Rashba spin precession and Datta-Das 
oscillations in spin-dependent transport are affected by e- 
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c interactions, Rashba SOI effects on electronic transport 
in interacting quantum wires have been studied in recent 
papers.— ^^^i 3 ^ 3 ^^" 3 - 7 - In effect, however, all those 
works only took e-e forward scattering processes into ac- 
count. Because of the Rashba SOI, one obtains a mod- 
ified LL phase with broken spin-charge separation ) 33 ' 34 
leading to a drastic influence on observables such as the 
spectral function or the tunneling density of states. Mo- 
roz et al. argued that e-e backscattering processes are 
irrelevant in the renormalization group (RG) sense, and 
hence can be omitted in a low-energy theory . 33 ' 34 Unfor- 
tunately, their theory relies on an incorrect spin assign- 
ment of the subbands ) 15 ' 22 which then invalidates several 
aspects of their treatment of interaction processes. 

The possibility that e-e backscattering processes be- 
come relevant (in the RG sense) in a Rashba quantum 
wire was raised in Ref. [38l . where a spin gap was found 
under a weak-coupling two-loop RG scheme. If valid, this 
result has important consequences for the physics of such 
systems, and would drive them into a spin-density-wave 
type state. To establish the spin gap, Ref. [H starts from 
a strict ID single-band model and assumes both a and 
the e-e interaction as weak coupling constants flowing 
under the RG. Our approach below is different in that 
we include the Rashba coupling a from the outset in the 
single-particle sector, i.e., in a nonperturbative manner. 
We then consider the one-loop RG flow of all possible in- 
teraction couplings allowed by momentum conservation 
(for not too small a) . This is an important difference to 
the scheme of Ref. |38|, since the Rashba SOI eliminates 
certain interaction processes which become momentum- 
nonconserving. This mechanism is captured by our ap- 
proach. The one-loop RG flow then turns out to be equiv- 
alent to a Kosterlitz-Thouless flow, and for the initial val- 
ues realized in this problem, e-e backscattering processes 
are always irrelevant. Our conclusion is therefore that no 
spin gap arises because of SOI, and a modified LL picture 
is always sufficient. We mention in passing that in the 
presence of a magnetic field (which we do not consider) , 
a spin gap can be present because of spin-nonconserving 
e-e "Cooper" scattering processesj 39 ' 40 the effects of e-e 
forward scattering in Rashba wires with magnetic field 
were studied as wel h 41 ' 42 ' 43 ' 44 Below, we also provide es- 
timates for the renormalized couplings entering the mod- 
ified LL theory, see Eq. (|26|) below. When taking bare 
(instead of renormalized) couplings, we recover previous 
results^ Note that the SOI in carbon nanotubes^ 5 - or 
graphenc ribbons 4 ^ leads to a similar yet different LL de- 
scription. In particular, for (achiral) carbon nanotubes, 
the leading SOI does not break spin-charge separation^ 5 . 
We here only discuss Rashba SOI effects in semiconduc- 
tor quantum wires in the absence of magnetic fields. 

We apply our formalism to a study of the Ruderman- 
Kittel-Kasuya-Yosida (RKKY) interactio n 47 ' 48 between 
two spin-1/2 magnetic impurities, £1,2, separated by 
a distance x. The RKKY interaction is mediated by 
the conduction electrons in the quantum wire which 
are exchange-coupled (with coupling J) to the impurity 



spins. In the absence of both the e-e interaction and 
the SOI, one finds an isotropic exchange (Heisenberg) 
Hamiltonian^ 8 - 



ffRKKY = - J 2 F ex (x) Si ■ £ 2 , F CK (x) (X 



cos(2fci?x) 



(2) 

where the 2fcir-oscillatory RKKY range function F ex (x) 
is specified for the ID case. When the spin SU(2) sym- 
metry is broken by the SOI, spin precession sets in and 
the RKKY interaction is generally of a more compli- 
cated (twisted) form. For a noninteracting Rashba quan- 
tum wire, it has indeed been establishe d 49 ' 50 ' 51 that the 
RKKY interaction becomes anisotropic and thus has a 
tensorial character. It can always be decomposed into 
an exchange (scalar) part, a Dzyaloshinsky-Moriya-likc 
(vector) interaction, and an Ising- like (traceless sym- 
metric tensor) coupling. On the other hand, in the pres- 
ence of e-e interactions but without SOI, the range func- 
tion has been showr* 5 ^ to exhibit a slow power-law decay, 
F e x(x) (x cos(2kFx)\x\~ ri , with an interaction- dependent 
exponent 77 < 1. The RKKY interaction in interacting 
quantum wires with SOI has not been studied before. 

For the benefit of the focussed reader, we briefly sum- 
marize the main results of our analyis. The effective low- 
energy theory of an interacting Rashba quantum wire is 
given in Eq. (|29[> . with the velocities (j30|) and the di- 
mcnsionless interaction parameters (pTTj) . Previous theo- 
ries did not fully account for the e-e backscattering pro- 
cesses, and the conspiracy of these processes with the 
broken SU(2) invariance due to spin-orbit effects leads 
to K s < 1 in Eq. (pTTj) . This in turn implies novel ef- 
fects in the RKKY interaction of an interacting Rashba 
wire. In particular, the power-law decay exponent in an 
interacting Rashba wire, see Eq. (f38|) . depends explic- 
itly on both the interaction strength and on the Rashba 
coupling. 

The structure of the remainder of this paper is as fol- 
lows. In Sec. HU we discuss the bandstructure. Inter- 
action processes and the one-loop RG scheme are dis- 
cussed in Sec. IIII1 while the LL description is provided 
in Sec. lIVl The RKKY interaction mediated by an inter- 
acting Rashba quantum wire is then studied in Sec. [V] 
Finally, we offer some conclusions in Sec. IVI1 Technical 
details can be found in the Appendix. Throughout the 
paper we use units where H = 1. 



II. SINGLE-PARTICLE DESCRIPTION 

We consider a quantum wire electrostatically confined 
in the z- direction within the 2DEG (xz-plane) by a har- 
monic potential, V c (z) = muj 2 z 2 /2, where m is the effec- 
tive mass. The noninteracting problem is then defined 
by the single-particle Hamiltonian 3 - ' 13 ' 14 ' 15 ' 17 



if si 



I 

2m 



(pi + pi) + V c {z) + a (a z p x - a x p z 



(3) 
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where a is the Rashba coupling and the Pauli matri- 
ces cr x z act in spin space. For a = 0, the transverse 
problem is diagonal in terms of the familiar ID har- 
monic oscillator eigenstates (Hcrmite functions) H n (z), 
with n = 0,1,2,... labeling the subbands (channels). 
Eigenstates of Eq. (J3j) have conserved longitudinal mo- 
mentum p x = k, and with the z-direction as spin quan- 
tization axis, cr z \cr) = a\a) with a =T,i= ±, the <r x p z 
term implies mixing of adjacent subbands with associ- 
ated spin flips. Retaining only the lowest (n = 0) sub- 
band from the outset thus excludes spin relaxation. We 
follow Ref. [HI and keep the two lowest bands, n = and 
n = 1. The higher subbands n > 2 yield only tiny cor- 
rections, which can in principle be included as in Ref. UtI 
The resulting 4x4 matrix representing H sp in this trun- 
cated Hilbert space is readily diagonalizcd and yields four 
energy bands. We choose the Fermi energy such that only 
the lower two bands, labeled by s = ±, are occupied, and 
arrive at a reduced two-band model, where the quantum 
number s = ± replaces the spin quantum number. The 
dispersion relation is 



2m 



UJ 

— + sak 



(4) 



with eigenfunctions ~ e lkx '4>k,s(z)- The resulting asym- 
metric energy bands ([4| are shown in Fig. [IJ The trans- 
verse spinors (in spin space) are given by 



t>k-(z) 



icos[0+(fc)]#i(3) 
sm[6+{k)}H {z) 

sm[9_(k)]H Q (z) 
icos[0_ (z) 



(5) 



with fc-dependent spin rotation angles (we take < 
s {k) < tt/2) 



.(ft) 



■ cot 



-2sk — to/a 
\f2mLU 



Os(-k). (6) 



As a result of subband mixing, the two spinor compo- 
nents of <j>k,s(z) carry a different z-dependence. They are 
therefore not just the result of a SU(2) rotation. For 
a = 0, we recover 9 S = tt/2, corresponding to the usual 
spin up and down eigenstates, with Hq(z) as transverse 
wavefunction; the s = + (s = — ) component then de- 
scribes the a ={ (a =f) spin eigenstate. However, for 
a^O, a peculiar implication of the Rashba SOI follows. 
From Eq. © we have \\m.k^±oo s (k) = (1 ± s)ir/4, such 
that both s = ± states have (approximately) spin a =1 
for k —> oo but a = j for k — > — oo; the product of spin and 
chirality thus always approaches asgn(fc) = — 1. More- 
over, under the time-reversal transformation, T = ia y C 
with the complex conjugation operator C, the two sub- 
bands are exchanged, 



-ikx , 



- k ,s{z) = sT[e^ M (z)], E_ s {-k) 



E s (k). 
(7) 
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FIG. 1: (Color online) Schematic band structure (j4j) of a typ- 
ical ID Rashba quantum wire. The red/blue curves show the 
s = ± bands, and the dotted curves indicate the next subband 
(the Fermi energy e.F is assumed below that band). For the 
low-energy description, we linearize the dispersion. It is no- 
tationally convenient to introduce bands A (solid lines) and 
B (dashed lines). Green and black arrows indicate the re- 
spective spin amplitudes (exaggerated). The resulting Fermi 



momenta are ±fc 



(A,B) 



with Fermi velocities va b ■ 



Time-reversal symmetry, preserved in the truncated de- 
scription, makes this two-band model of a Rashba quan- 
tum wire qualitatively different from Zecman-spin-split 
modcls<2I 

In the next step, since we are interested in the 
low-energy physics, we linearize the dispersion relation 

around the Fermi points ±k^' B \ see Fig.[TJ which results 
in two velocities va and vg, see Eq. ([TJ) . The lineariza- 
tion of the dispersion relation of multi-band quantum 
wires around the Fermi level is known to be an excellent 
approximation for weak e-e interactions^ Explicit val- 
ues for S in Eq. ((T|) can be derived from Eq. (j4|), and we 
find 5(a) cx a 4 for a — > 0, in accordance with previous 
estimates^ We mention that 6 < 0.1 has been estimated 
for typical geometries in Ref. [3J. The transverse spinors 
<j>k 8 {z), Eq. (0, entering the low-energy description can 

be taken at k = ±kp B \ where the spin rotation angle 
© only assumes one of the two values 



(8) 



The electron field operator ^(x, z) for the linearized two- 
band model with v = A, B = +, — can then be expressed 
in terms of ID fermionic field operators ip v ^ r (x), where 
r = R,L = +,— labels right- and left-movers, 



E 1 

v.r— ± 



(9) 



with cf>k,s{z) specified in Eq. ([5]). Note that in the left- 
moving sector, band indices have been interchanged ac- 
cording to the labeling in Fig. [T] 

In this way, the nonintcracting second-quantized 
Hamiltonian takes the standard form for two incquiva- 
lcnt species of ID masslcss Dirac fcrmions with different 



4 



velocities, 



Hq = -i ^2 rv v I dx ipt ir d x i/) Vt , 



(10) 



The velocity difference implies the breaking of the spin 
SU(2) symmetry, a direct consequence of SOI. For a = 0, 
the index v coincides with the spin quantum number a 
for left-movers and with —a for right-movers, and the 
above formulation reduces to the usual Hamiltonian for 
a spinful single-channel quantum wire. 



condition (|13[) eliminates one important interaction pro- 
cess available for a = 0, namely interband backscattering 
(see below) . This is a distinct SOI effect besides the bro- 
ken spin SU{2) invariance. Obtaining the complete "g- 
ology" classification^ of all possible interaction processes 
allowed for a ^ is then a straightforward exercise. The 
corresponding values of the interaction matrix elements 
are generally difficult to evaluate explicitly, but in the 
most important case of a thin wire, 



1 



(16) 



III. INTERACTION EFFECTS 

Let us now include e-e interactions in such a single- 
channel disorder-free Rashba quantum wire. With the 
expansion @ and r = (x,z), the second-quantized two- 
body Hamiltonian 

Hi = lj *t( ri )*t( ra )y( ri _ r 2 )*(r 2 )*(n) 

(11) 

leads to ID interaction processes. We here assume 
that the e-e interaction potential V(yi — r 2 ) is exter- 
nally screened, allowing to describe the ID interactions 
as effectively local. Following standard arguments, for 
weak e-e interactions, going beyond this approximation 
at most leads to irrelevant corrections.— We then obtain 
the local ID interaction Hamiltonian— 

Hl = \Yl V i^n} J dx ^i,r 1 V4,r 2 VV 3 ,r 3 ' ! /W4> 

(12) 

where the summation runs over all quantum numbers 
v\, . . . ,va and r\, . . . , subject to momentum conserva- 
tion, 



nk { p ] + r 2 k { p ] = r 3 k ( p ] + r 4 fc 



F ■ 



(13) 



With the momentum transfer q = r-ykp — r&kp and 
the partial Fourier transform 



V(q;z) = / dx e~ iqx V{x,z) 



(14) 



of the interaction potential, the interaction matrix ele- 
ments in Eq. (|12|) are given by 



dzidz 2 V(q; Z\ - z 2 ) 

At 



j>4) 



(Zl) 

(**)• (15) 



Since the Rashba SOI produces a splitting of the Fermi 
momenta for the two bands, k^ — kp ~ 2am, the 



where d is the screening length (representing, e.g., the 
distance to a backgate), analytical expressions can be 
obtained^ To simplify the analysis and allow for ana- 
lytical progress, we therefore employ the thin-wire ap- 
proximation (|16p in what follows. In that case, we can 
neglect the z dependence in Eq. (fl"4"|) . Going beyond this 
approximation would only imply slightly modified values 
for the e-e interaction couplings used below. Using the 
identity 



dz 

8„,ji 8. 



(z) = 



(17) 



where the angles 9a, b were specified in Eq. only two 
different values Wq and W\ for the matrix elements in 
Eq. I|15p emerge. These nonzero matrix elements are 



V ur ,v'r',v>r',vr = W = V(q = 0), 
'i/r,i/ , r' ) — i/'— r', — v— r = *v\ 



(18) 



— 6b) V (q = kp 



k 



(B) 



We then introduce ID chiral fcrmion densities p vr {x) = 
: iptr^fr '■■> where the colons indicate normal-ordering. 
The interacting ID Hamiltonian is H = Hq + Hj with 
Eq. (JTOD and 

Hi = ^ / ' dx([g 2 \\ v 5„y + g2±8 v ,-u']8 r - r ' 

vv' ,rr' 

+ \94,\\v8v,V + g4J.8 v ,-u']8 r ,r'J PvrPu'r 1 (19) 

The e-e interaction couplings arc denoted in analogy to 
the standard g-ology, whereby the (74 ((72) processes de- 
scribe forward scattering of ID fermions with equal (op- 
posite) chirality r = R, L = +,—, and the labels ||, _L, 
and / denote intraband, interband, and band flip pro- 
cesses, respectively. Since the bands v = A, B = +, — are 
inequivalent, we keep track of the band index in the intra- 
band couplings. The gf term corresponds to intraband 
backscattering with band flip. The interband backscat- 
tering without band flip is strongly suppressed since it 



5 



does not conserve total momentum^ and is neglected in 
the following. For a = 0, the fl'.o/j, couplings coincide 
with the usual ones^ 2 - for spinful electrons, while <?/ re- 
duces to <7ij_ and 92,\\/± — * 92,±/\\ due to our exchange 
of band indices in the left-moving sector. According to 
Eq. (fl8|) . the bare values of these coupling constants are 



9a\\ v = 941- = 9i\\ v = W , 
92±=W -Wi, gf = Wi. 



(20) 



The equality of the intraband coupling constants for the 
two bands is a consequence of the thin-wire approxima- 
tion, which also eliminates certain exchange matrix ele- 
ments. 

The Hamiltonian Ho + Hi then corresponds to a spe- 
cific realization of a general asymmetric two band-model, 
where the one-loop RG equations are know m 54 ' 57 Using 
RG invariants, we arrive after some algebra at the two- 
dimensional Kosterlitz-Thoulcss RG flow equations, 



dm 

dl 



-2 d 9f 



(21) 



for the rescaled couplings 

9i\\A 9i\\ 



!)2 



9f 



2itva 2itvb kvf 



(22) 



1 + 7 3/ 



2 ttvf 

where we use the dimensionless constant 



1 



7 



V A VB 



1-S 2 



> 1. 



(23) 



As usual, the 54 couplings do not contribute to the one- 
loop RG equations. The initial values of the couplings 
can be read off from Eq. (f20|) , 



m = 0) 
9f(l = 0) 



(7 - 1)W + Wi 
ttvf 



1 + 7 Wi 

2 TTVp 



(24) 



The solution of Eq. (|2"Tj) is textbook material,— and 5/ 
is known to be marginally irrelevant for all initial condi- 
tions with |.g/(0)| < ff 2 (0). Using Eqs. |[T8] ) and (|24 j) . this 
implies with 7 ~ 1 + S 2 the condition 



V(0) > -cos 2 (^ 



) V ( k ( F A) + k 



SB) 



(25) 



which is satisfied for all physically relevant repulsive e- 
e interaction potentials. As a consequence, intraband 
backscattering processes with band flip, described by the 
coupling g~f, are always marginally irrelevant, i.e., they 
flow to zero coupling as the energy scale is reduced, 
9f = 9f{l ~ * °°) = 0. Therefore no gap arises, and a 
modified LL model is the appropriate low-energy theory. 
We mention in passing that even if we neglect the ve- 
locity difference in Eq. JT]), no spin gap is expected in 



a Rashba wire, i.e., the broken SU{2) invariance in our 
model is not required to establish the absence of a gap. 

The above RG procedure also allows us to extract 
renormalized couplings entering the low-energy LL de- 
scription. The fixed-point value g\ = ,92 (I — * 00) now 
depends on the Rashba SOI through 7 in Eq. (J23J). With 
the interaction matrix elements Wo,i m Eq. (fTS)) . it is 
given by 



v /[(7-l)^o + lU 1 ] 2 -(7 + l)^f72 
92 = ■ (26) 



TTVF 



For a = 0, we have 7 = 1 and therefore g\ = 0. The 
Rashba SOI produces the nonzero fixed-point value (|26[) . 
reflecting the broken SU (2) symmetry. 



IV. LUTTINGER LIQUID DESCRIPTION 

In this section, we describe the resulting effective 
low-energy Luttinger liquid (LL) theory of an inter- 
acting single-channel Rashba wire. Employing Abclian 
bosonization^ we introduce a boson field and its con- 
jugate momentum for each band v = A, B = +,—. 
It is useful to switch to symmetric ("charge"), $ c (cc) 
and H c (x) — —d x O c (x), and antisymmetric ("spin" for 
a = 0), <fr s (x) and II s (a;) = — d x Q s , linear combinations 
of these fields and their momenta. The dual fields $ 
and O then allow to express the electron operator from 
Eq. (JHJ) and the "bosonization dictionary," 



x e 



rk\. .vr 



c+i v '7r/2[r$ c -|-e c +i/r* 3 +i/e s ] 



(27) 



where a is a small cutoff length and the stan- 

dard Klein factorsi 32 i 52 ' 58 (To recover the conventional 
expression for a = 0, due to our convention for the 
band indices in the left-moving sector, one should re- 
place $ S ,6 S — ► -6 S! -$ S .) Using the identity fj"7|). we 
can now express the ID charge and spin densities, 



p{x) = J dz&V, S(x) = J dztf^V, 



(28) 



in bosonized form. The (somewhat lengthy) result can 
be found in Appendix [A] 

The low-energy Hamiltonian is then taken with the 
fixed-point values for the interaction constants, i.e., 
backscattering processes are disregarded and only appear 
via the renormalized value of g\ in Eq. (|26[) . Following 
standard steps, the kinetic term Hq and the forward scat- 
tering processes then lead to the exactly solvable Gaus- 
sian field theory of a modified (extended) Luttinger liq- 
uid, 



H = E v iJ dx U^ + ^id^A (29) 

j—c,s 3 s 

+ v x J dx (^k x u c u s + -L(d x $ c )(d x <s> s ^j . 
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Using the notations 54 



Vs 

in 



- Wq/ttvf and 

A ~ 92\\B 



92\\A + 9* 2 \\B ± 2 9*2± 



where explicit (but lengthy) expressions for the fixed- 
point values g2\\A/B an( i 92± can be straightforwardly 
obtained from Eqs. (f!22")) and (f2"u| . the renormalized ve- 
locities appearing in Eq. (|29|) are 



v F \J{^ + 9a) 2 -y 



W \ 2 [2W a -W^ 2 




ttvf J V 2ttvf 
• v F , 



(30) 



In the respective second equalities, we have specified the 
leading terms in \6\ <C 1, since the SOI- induced relative 
velocity asymmetry 5 is small even for rather large a, 
see Eq. (JTJ. The corrections to the quoted expressions 
are of 0(S 2 ) and are negligible in practice. It is note- 
worthy that the "spin" velocity v s is not renormalized 
for a Rashba wire, although it is well-known that v s will 
be renormalized due to W\ for a = This difference 
can be traced to our thin- wire approximation Q16p . When 
releasing this approximation, there will be a renormaliza- 
tion in general. Finally, the dimensionless LL interaction 
parameters in Eq. ([25)) are given by 



K r = 



K s = 



K x = 



.94 + y+ 




2nv F + Wi 
2irv F + 4 Wo - Wi 



1 \6\, 



V2 



TTVF 



(31) 



A-kvf + W\ 
Attvf — Wi 



where the second equalities again hold up to contribu- 
tions of 0(S 2 ). When the 2k f component of the in- 
teraction potential W\ = 0, see Eq. (Tig)) , we obtain 
K s = K\ = 1, and thus recover the theory of Ref. [22|. 
The broken spin SU (2) symmetry is reflected in K s < 1 
when both (5^0 and W\ 7^ 0. 

Since we arrived at a Gaussian field theory, Eq. ([25]) . 
all low-energy correlation functions can now be computed 
analytically without further approximation. The linear 
algebra problem needed for this diagonalization is dis- 
cussed in App. [A] 



V. RKKY INTERACTION 

Following our discussion in Sec. HI we now investigate 
the combined effects of the Rashba SOI and the e-e in- 
teraction on the RKKY range function. We include the 
exchange coupling, H' = JJ2i=i 2 ^« ' S( x i)> °f tne ID 
conduction electron spin density S(x) to localized spin- 
1/2 magnetic impurities, separated by x = x\ — x%- The 
RKKY interaction Hrkky, describing spin-spin interac- 
tions between the two magnetic impurities, is then ob- 
tained by perturbation theory in Ji 4 ^ In the simplest ID 
case (no SOI, no interactions), it is given by Eq. ©. In 
the general case, one can always express it in the form 



RKKY 



-J 2 ]TF Qb (x)£^, 



(32) 



with the range function now appearing as a tensor (/3 
1/ksT for temperature T), 



F a \x) 



dT X a \x,T). 



(33) 



Here, the imaginary-time (r) spin-spin correlation func- 
tion appears, 



X ab (x,r) = (S a (x,T)S b (0,0)). 



(34) 



The ID spin densities S a (x) (with a = x,y,z) were 
defined in Eq. Q28|) . an( l their bosonized expression is 
given in App. [Al which then allows to compute the cor- 
relation functions (f34|) using the unperturbed (J = 0) 
LL model (f29|) . The range function thus effectively co- 
incides with the static space-dependent spin suscepti- 
bility tensor. When spin SU(2) symmetry is realized, 
X ab (x) = 8 ab F ex (x), and one recovers Eq. ([2]), but in gen- 
eral this tensor is not diagonal. For a LL without Rashba 
SOI, F ex (x) is as in Eq. @ but with a slow power-law 
decay^ 

If spin SU(2) symmetry is broken, general arguments 
imply that Eq. (|32|) can be decomposed into three terms, 
namely (i) an isotropic exchange scalar coupling, (ii) a 
Dzyaloshinsky-Moriya (DM) vector term, and (iii) an 
Ising-like interaction, 

#rkky/J 2 = -F ex (x)-£i • £ 2 - F BM (x) ■ (£1 X S 2 ) 

(35) 



a . b 



where F cx (x) = i^^W- The DM vector has the 
components 



and the Ising-like tensor 




cccab 



3 ^ 



F cc 8 



(x) 
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is symmetric and tracclcss. For a ID nonintcracting 
quantum wire with Rashba SOI, the "twisted" RKKY 
Hamiltonian |35|) has recently been discussed , 49 ! 50 ' 51 and 
all range functions appearing in Eq. Q35p were shown to 
decay oc |a;| _1 , as expected for a noninteracting system. 
Moreover, it has been emphasized 5 - - that there are dif- 
ferent spatial oscillation periods, reflecting the presence 
of different Fermi momenta k^' B ^ in a Rashba quantum 
wire. 



Let us then consider the extended LL model (|29l) . 
which includes the effects of both the e-e interaction 
and the Rashba SOI. The correlation functions obey 
X ba {x,r) = x ab (— x, — t), and since we find \ xz = X yz 
. the anisotropy acts only in the xy-plane. The four 
nonzero correlators are specified in App. El where only 
the long-ranged 2kp oscillatory terms are kept. These 
are the relevant correlations determining the RKKY in- 
teraction in the interacting quantum wire. We note that 
in the noninteracting case, there is also a "slow" oscilla- 
tory component, corresponding to a contribution to the 



k {A) - k 

F "'F 



(B] 



Re- 



RKKY range function oc cos 

markably, we find that this 1 jx decay law is not changed 
by interactions. However, we will show below that inter- 
actions cause a slower decay of certain "fast" oscillatory 
terms, e.g., the contribution oc cos(2fc] ? x). We there- 
fore do not further discuss the "slow" oscillatory terms 
in what follows. 

Collecting everything, we find the various range func- 
tions in Eq. (|35|) for the interacting case, 



F ex (x) = 
+ 

Fdm(x) = 
(x) = 



i [(1 + cos 2 (2^)) cos (2kPx) fW(x) 

V 

cos 2 {9 A + 9 B )cos \(k ( p ] + k ( p ] )x] F^\x)] , 
e^^cos(2^)sin hk^x) F^(x), (36) 



iY,Gl(x)-F ex (a 



with the auxiliary vector 



/ 



cos" 



cos (2ft£° a:) (x) \ 
cos 2 (2^)cos (2kPx) F^\x) 
9 A + B ) cos + k { p } )x\ f!) 2) (x) , 



(1 2") 

The functions F„ ' (x) follow by integration over r from 
F^' 2) (x,t), see Eqs. (jATJ and JA2|) in App. El This 
implies the respective decay laws for a <gc |x| <gC vp/ksT, 

(37) 



(x) oc |a/a;|" 
Fj?Xx) oc \a/x\- 1+K ' +1 ' K ' 



All those exponents approach unity in the noninteracting 
limit, in accordance with previous results! 49 ' 50 Moreover, 



in the absence of SOI [a. = 6 = 0), Eq. pTj) reproduces 
the known |x| _is:c decay law for the RKKY interaction 
in a conventional LL^ 2 - 

Since K s < 1 for an interacting Rashba wire with 
5^0, see Eq. (|3"Tj) . we conclude that F^ with v = B, 
corresponding to the slower velocity vb = vf(1 — S), leads 
to the slowest decay of the RKKY interaction. For large 
distance x, the RKKY interaction is therefore dominated 

( B) 

by the 2k F oscillatory part, and all range functions de- 
cay oc |x|~ ,,B with the exponent 



f] B =K C + K„-1-2[1 



K, 



K{) v c + v. 



< 1. (38) 



This exponent depends both on the e-e interaction po- 
tential and on the Rashba coupling a. The latter de- 
pendence also implies that electric fields are able to 
change the power-law decay of the RKKY interaction in 
a Rashba wire. The DM vector coupling also illustrates 
that the SOI is able to effectively induce off-diagonal cou- 
plings in spin space, reminiscent of spin precession ef- 

( B) 

fects. Also these RKKY couplings are 2k F ' oscillatory 
and show a power-law decay with the exponent 



VI. DISCUSSION 

In this paper, we have presented a careful derivation of 
the low-energy Hamiltonian of a homogeneous ID quan- 
tum wire with not too weak Rashba spin-orbit inter- 
actions. We have studied the simplest case (no mag- 
netic field, no disorder, single-channel limit), and in par- 
ticular analyzed the possibility for a spin gap to oc- 
cur because of electron-electron backscattering processes. 
The initial values for the coupling constants entering the 
one-loop RG equations were determined, and for rather 
general conditions, they are such that backscattering is 
marginally irrelevant and no spin gap opens. The re- 
sulting low-energy theory is a modified Luttingcr liquid, 
Eq. (|2T)]) . which is a Gaussian field theory formulated 
in terms of the boson fields <& c (x) and Q s (x) (and their 
dual fields). In this state, spin-charge separation is vi- 
olated due to the Rashba coupling, but the theory still 
admits exact results for essentially all low-energy corre- 
lation functions. 

Based on our bosonized expressions for the ID charge 
and spin density, the frequency dependence of various 
susceptibilities of interest, e.g., charge- or spin-density 
wave correlations, can then be computed. As the calcu- 
lation closely mirrors the one in Refs. I3411351 . we do not 
repeat it here. One can then infer a "phase diagram" 
from the study of the dominant susceptibilities. Accord- 
ing to our calculations, due to a conspiracy of the Rashba 
SOI and the e-e interaction, spin-density-wave correla- 
tions in the xy plane are always dominant for repulsive 
interactions. 

We have studied the RKKY interaction between two 
magnetic impurities in such an interacting ID Rashba 
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quantum wire. On general grounds, the RKKY interac- 
tion can be decomposed into an exchange term, a DM 
vector term, and a traceless symmetric tensor interac- 
tion. For a noninteracting wire, the corresponding three 
range functions have several spatial oscillation periods 
with a common overall decay oc | a; | 1 . We have shown 
that interactions modify this picture. The dominant con- 
tribution (characterized by the slowest power-law decay) 

( B) 

to the RKKY range function is now 2k F oscillatory for 
all three terms, with the same exponent t\b < 1, see 
Eq. ([55)) . This exponent depends both on the interac- 
tion strength and on the Rashba coupling. This raises 
the intriguing possibility to tune the power-law exponent 
rjB governing the RKKY interaction by an electric field, 
since a is tunable via a backgate voltage. We stress again 
that interactions imply that a single spatial oscillation 

( B) 

period (wavelength n/k F ) becomes dominant, in con- 
trast to the noninteracting situation where several com- 
peting wavelengths are expected. 

The above formulation also holds promise for future 
calculations of spin transport in the presence of both in- 
teractions and Rashba spin-orbit couplings, and possibly 
with disorder. Under a perturbative treatment of impu- 
rity backscattering, otherwise exact statements are pos- 
sible even out of equilibrium. We hope that our work will 
motivate further studies along this line. 
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APPENDIX A: BOSONIZATION FOR THE 
EXTENDED LUTTINGER LIQUID 

In this appendix, we provide some technical details re- 
lated to the evaluation of the spin-spin correlation func- 
tion under the extended Luttinger theory (f2"5|) . The ex- 
act calculation of such correlations is possible within the 
bosonization framework, and requires a diagonalization 
of Eq. (j29l). 

The ID charge and spin densities (|2"8|) can be written 
as the sum of slow and fast (oscillatory) contributions. 
Using Eq. (fT?)) . the bosonized form for the ID charge 
density is 



p(x) 



2 * ^ 2l t 

- Ox®c VARVAL COS( 

7r ira 



x sin 



2tt$ c 



Similarly, using the identity 



dz 



cos(v27r9 s ) 



cos (0a - 6 B )S v -u' 

S r ,r' ( — IVT COS {6 A ~\~ @b) ^v,—v' 

vr cos (2 

+5 r - r i \ —ivr cos (2 6 V ) 5 v y 
vr cos (6a + Ob) 5 v - v > 



the ID spin density vector has the components 



ox/ \ -VarVbr , n 

S (x) = —% cos (0a 

na 



■ COS 



/ B )cos[(4 A) -4 B) )^ + 



2tt$ s 



27T$ r 



Sill 



sin(V27r6 s ) 



C y i \ -VARVBR , a . \ . 
S y (x) — l COS (0A + 0_b) Sill 

ira 

-« £ 

V= A.B=+,- 



J,(A) ,(By 
F F ' 



2tt$ s 



sin(v / 2~7r6 s 



v cos(20„) cos 

2ita 



2tt ($ c + z/$ s 



S'(x) 



1 



V8tt 
.VarVbl 



[(cos 20 A + cos 20 B ) d x e s + (cos 20 A - cos 20 B ) d x Q c ] 



COS(0A + 0b) COS 



1A) AB) 



27T$ f 



sin(V27r$ s ) 



Note that while d x $ c is proportional to the (slow part of by both c and s sectors, 
the) charge density, the (slow) spin density is determined 

Next we specify the nonzero components of the 
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imaginary-time spin-spin correlation function x ab ( x i T )-, an d 
see Eq. (fMf . Using the above bosonized expressions, 
some algebra yields 

cos ( 2kp S> x ) 



iscos(29v) sin (^k^x^J 

cos 2 (2^)cos( 2 A>U . XXV(X,T) = 5 2(2^ ^W)- 



x»(*,r)=i; ^ — t-m*,r), 



X zz (x lT ) = 



^-v COS (W^i 



1 1 fcp + & 



2(27ra) 2 



[(» 



~ (1 2) 

Here, the functions F^ A B=+ 



(x, t) are given by 



h i hx,t)= n 

3=1,2 



fiuj . / n(ujT — ix) 



sin 



and 



^ {x,r)= n 

3=1,2 



/3iij 



_( T U) +T U) 



I Tr(ujT-ix) \ 



sin 



The dimensionless numbers T^> appearing in the expo- 
nents follow from the straightforward (but lengthy) diag- 
onalization of the extended LL Hamiltonian (|29|) . where 
the Uj are the velocities of the corresponding normal 
modes. With the velocities ([50)) and the dimensionless 
Luttinger parameters (piTj) , the result of this linear alge- 
bra problem can be written as follows. The normal-mode 
velocities u\ and ui are 



2U 3 2 =1,2 



Kl 
KJC t 



KcKs 















)+vl- 







1/2 



and the exponents 1 ' 2 ) appearing in F„ (x, r) arc 
given by 



-.0') 



-.0') 



Uj(uf-ti|) 
(-l)^ A f A 



J K c K s v c 



J K c K s v s 



K C K S 



e s e 3 



(-l) j Vs 



r 



K s Uj(ul ~ u i) 



K c K s v\v c ' 



0) 



A s 



a; 



Since | (5 1 <C 1, we now employ the simplified expres- 
sions for the velocities in Eq. (f50|) and the Luttinger liq- 
uid parameters in Eq. (f3Tj) . which are valid up to 0(S 2 ) 
corrections. In the interacting case, this yields for the 
normal- mode velocities simply u\ = v c and 1*2 = v s . 
(In the noninteracting limit, the above equation instead 
yields Ui = va and U2 = t's, see Eq. ([1]).) Moreover, the 
exponents simplify to 



r 



(i) 



Kr, r 



(2) 



= r 



(i) 



_ r d) _ 



p(2) _ ^ p(2) 



,(2) 



''A 



J C s 



(K x v c + K c v s /K x ), 
:(K x v s + K c v c /K x ), 
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F (1.2) _ , r (2) 

1 $ c e s — 11=1 $ c $ 3 



and 



Collecting everything and taking the zero-temperature 

~ (1 2) 

limit, the functions F„J.± (x, r) take the form 



V C T — IX 



W S T — IX 



(Al) 









V C T — ix 


V S T — ix 


a 




a 



-1/K, 



v x (K x v s + K c v c /K x ) 
~ V v'i-,^ 



(v s t — ix)(v c T + ix) 
{v s t + ix)(v c T — ix) 



The known form of the spin-spin correlations in a LL 
with a = is recovered by putting v\ oc S = 0. 
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